Mesoscopic Charge Density Wave in a Magnetic Flux 
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The stability of a Charge Density Wave (CDW) in a one-dimensional ring pierced by a Aharonov- 
Bohm flux is studied in a mean-field picture. It is found that the stability depends on the parity of 
the number N of electrons. When the size of the ring becomes as small as the coherence length f, 
the CDW gap increases for even N and decreases for odd N. Then when TV is even, the CDW gap 
decreases with flux but it increases when TV is odd. The variation of the BCS ratio with size and 
flux is also calculated. We derive the harmonics expansion of the persistent current in a presence 
of a finite gap. 
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I. INTRODUCTION 

It has been proposed recently that a Aharonov-Bohm 
(AB) flux <j> should affect the stability of a Charge Density 
Wave (CDW) in a one-dimensional ring geometrymrtJ 
it may have been seen in a recent experiment on CD 
pierced by AB flux lines trapped in columnar defectsl 
One possible cause for the modulation of the CDW sta- 
bility is the discreteness of the spectrum and the moduL 
tion of the position of the energy levels with the fluxH' 
Indeed, it has been predicted that the CDW gap and 
the critical temperature oscillate with the flux <f>, with 
a period cf>o — h/e and that they are maxima at <p = 
and minima at 4> = 4>a/2. When the perimeter L of the 
ring becomes of the order of the correlation length £, i.e. 
when the CDW order parameter becomes of the order of 
the mean level spacing, the CDW can even be.-destroyed 
when \(j> — n(f>o\ is larger than a critical valueclo. 

In this paper, we elaborate on these ideas and we show 
that the stability of the CDW depends crucially on the 
parity of the number N of particles in the ID ring (con- 
sidering here spinless particles), an effect which has not 
been properly considered in previous works. When this 
number N is odd, the CDW is destabilized by decreas- 
ing the length of the ring. But when N is even, it we 
find that the CDW is stabilized. When N is even, the 
effect oi an AB flux is to destabilize the CDW, as found 
in refs.0. But when N is odd, we find that the CDW 
can he^siabilized by the flux, contrary to the conclusions 
of refsJafl. This effect is reminiscent of the parity 
for the persistent current in one-dimensional ring: 
In the following, we shall consider the case of spinless 
electrons. 

In the next section, we establish the thermodynamic 
equations for the CDW in a finite ID system, with em- 
phasis on the parity effect. In section III, we calculate 
the flux variation of the order parameter in the form of 
a Fourier series. In section IV, the same is done for the 
critical temperature. Finally we calculate the persistent 
current in the last section. 
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II. CDW IN A MAGNETIC FLUX 

Consider a one-dimensional ring of perimeter L. The 
periodic boundary conditions fix the wave vector of the 
eigenstates. In the presence of a flux, the possible wave 
vectors are: 

k = — (p + ip) 

where ip is the dimensionless flux 4>/(f>o and p £ Z. Let N 
be the total number of electrons. The Fermi wave vector 
&f is: 

Nir 
k F = — 

and the Fermi energy e F is: 



e F = 



h 2 /A^7r x 2 
2m I ~~L 



(1) 



It is important to stress that with this choice of k F and 
tp, the number of electrons is independent of the flux, 
Fig.®- 
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FIG. 1. Dispersion relation in the metallic phase, a) when 
the number N of particles is even (do not forget the state 
k — 0), b) when TV is odd. The horizontal line gives the 
position of the Fermi level. 

Although the chemical potential is in principle fixed 
by a reservoir in ref.EI, we shall see later that their results 
correspond actually to a fixed even number of particles. 

The theory of the CDW ordering is well known. In the 
presence of a periodic potential with wave vector Q, the 
Hamiltonian is: 

H = Ho + 2A cos(Qr + 0) 

where the eigenenergies of Ho are known. Because of 
the periodic potential, the states k are in principle cou- 
pled to states k ± nQ where n = 1, 2, 3, We consider 

the case of a weak potential. Neglecting commensura- 
bility effects, one has to diagonalize a 2 x 2 matrix: the 
states with wave vector k ~ k F are coupled to states with 
wave vector k — Q ~ — kp\ 



e(k) Ae i6 
Ae- ie e(fc-Q) 



The eigenvalues are given by E — €fc+ ^ fc ± 



level spacing becomes as large as the the CDW gap, see 
Fig.(l). 



Ch e 2 k ° j +A 2 . It is then convenient to linearize 

the dispersion near the Fermi level so that, by taking the 
origin of the energies at the Fermi level, one has: 



TV 

e k = hv F (\k\ - hp) = 26 [ \p + ip\ - — 



where the mean level spacing 6 near the Fermi level is 
given by irhvp/L. Writing the nesting vector Q in the 
form Q — 2nq/L, where q is an integer, the eigenvalues 
in the CDW phase are: 




FIG. 2. Dispersion relation near the Fermi level, a) when 
the number iV of particles is even , b) when TV is odd. The 
horizontal line gives the position of the Fermi level. The small 
dots denote the states in the metallic phase and the large dots 
represent the states in the CDW phase 



E p = {q-N)5±\J^(p-^ + V y + ^ . 

The nesting condition implies q = N so that the nesting 
vector is 



Q = 2k F 



2Nn 
L 



and the energy levels in the CDW phase become: 



E, 



N 



,-±^45^-- 



tp) 2 + A 2 



N 

= ±25\j(p--+p) 2 + K 2 



(2) 



where A = A/(2<5). A measures the gap amplitude in 
units of the mean level spacing and can be written as 
the ratio of the length over the CDW coherence length 
£ = hvp/irA. One has A = L/(2tt 2 £). The effects dis- 
cussed in this paper appear for small rings when the mean 



From minimization of the total free energy 
A 2 

F = —L - k B TLogZ 
A 

where A is the interaction parameter and Z is the 
grand canonical partition function, the well-known self- 
consistency condition is obtained: 



where g = \/2Trhvp. When N is even, the introduction 
of a new variable n = p — N/2 reduces this equation to: 



1 = 5 



^tanh(4V(n + ^) 2 +A 2 ) 



2 N /(^ + y) 2 +A 2 



(3) 



where f3' — 2(35. This is the result found in ref.H. How- 
ever, when the number N is odd, n = p — (N — l)/2 and 
the self-consistency equation becomes: 



2 



1=3 



tanh(£-^/ (n 



<p) 2 + A 2 ) 



2 v /(n-i + (^) 2 +A 2 



(4) 



Thus the stability of the CDW depends on the parity of 
the number of electrons. The flux dependence for the 
even and odd parities are deduced from each other by a 
translation of </>o/2. We stress the fact that the results of 
ref.13 correspond to a fixed even number of particles, once 
the chemical potential is fixed to the value of eq. ([!]). If 
the chemical potential was fixed to any other value, the 
number of particles would vary with the field leading to 
discontinuities of the persistent currents and of the CDW 
stability. 



III. GROUND STATE 

We first study the evolution of the CDW order param- 
eter A v at zero temperature, as a function of the size and 
of the AB flux. When A~ is even, it is given by 



1 



E 



i 



^(n + ^) 2 + A 2 v 



(5) 



where the dimensionless parameter A v = A V /2S has 
been introduced. This sum diverges at large n and is 
usually cut-off at an energy scale E* of the order of the 
bandwidth. This corresponds to an integer n* = E* /26. 
In ref.EI, the flux dependence of the order parameter is 
given by: 



J{n + if) 2 + A% y/n 2 + A 2 







(G) 



which corresponds to the case where N is even. However, 
when N is odd, n should be replaced by n — 1/2 as seen 
in eq.(|). 

We found it convenient to write the Fourier decompo- 
sition of the flux dependence of the order parameter. To 
transform eq. (^J) , we use the Poisson summation formula 

^ f(n + <p) = go + 2 ^ 9m cos{2nmip) 



m>0 



where 



/>oo 

9m = 2 / f(y) cos(2nmy)dy . 
Jo 

The constant term go is divergent and must be cut-off 
with n* . The Fourier components are convergent. As a 
result one gets: 

- = In %— + 2 V (-l) Arm AT (27rmA ( ,) cos27rmv? (7) 
y v m>0 



The factor (— \} Nm in the harmonics expansion takes into 
account the parity |effcct. A'o is a modified Bessel function 
of the second kindB. 

It is instructive to consider first the ring in the absence 
of external flux. For the infinite bulk system, the gap A& 
is given by: 



1 _ 2E* 

~g ~ n 



(8) 



When the size becomes finite, the gap in zero flux Ao is 
given by: 

ln X £ = 2 E(- 1 ) Ar "^o(2^A b ) (9) 

b m>0 b 

where the parameter A^ = At/26 measures the bulk gap 
in units of mean level spacing and is proportional to the 
size. 

The variation of the gap with Af, = L /2-k 2 ^ is shown 
on Fig.(|J). £{, is the coherence length of the infinite sys- 
tem. When the number of particles is even, one sees that 
the gap increases with decreasing size. This seems to be 
in contradiction with reffl who found a decrease of the 
critical temperature T c when the size decreases. We will 
comment on this point in the next section. When the 
number of particles is odd, the gap decreases with de- 
creasing size. When the length becomes too small, the 
order parameter can even vanish. This happens for a crit- 
ical value of Ah = I/27 = .280. 7 is the Euler constant 
7 = 1.781. 



2.5 



Ao 



1 

0.5 





0.2 



0.4 



0.0 



0.8 



A 



FIG. 3. Variation of the dimensionless gap Ao/Aj, with 
At — At/28. At is the gap of the infinite system. A;, is 
proportional to the size L. 

We now turn to the effect of the AB flux. In a finite 
flux, the gap equation (^J) transforms into: 

ln^ = 2 ^(-l) JVm Ar (27rm^A fc )cos27rmv 3 (10) 

k m>0 k 

Fig.(^) shows the gap in the case of even N, for different 
fluxes. It is seen the effect of the flux is to reduce the gap, 
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this effect being larger for small sizes. Actually the case 
ip = 1/2 with an even number of particles is equivalent 
to the case (p = with an odd number of particles. This 
is obvious in Figs.(^) and (^) and in the structure of the 
gap equation (j?]). 
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FIG. 4. Variation of the dimensionless gap A v / A^ with At, 
for different fluxes tp = 0, 0.1, 0.2, 0.3, 0.4, 0.5 . 

Fig.(||) shows the flux dependence of A^/Aj, with the 
flux for different parameters Ab/26 cx L/^. For an even 
number of particles, the order parameter is enhanced at 
zero flux and is reduced for large flux. When Ab/28 
becomes too small, the order parameter can even van- 
ish near tp = 1/2. This happens for a critical value of 
A;, = I/27 = .280. For smaller rings, i.e. when A;, is 
smaller, the CDW disappears at a critical flux given by: 

iP{ip c )+ip(l-<p c ) =21nA 6 /2 

where ip is the digamma function. 



1.4 
1.2 
1 

A 9 0.8 
A b 0.6 
0.4 
0.2 











•5 






.35 . 


Even 








A=.25\ 

b \ 





0.1 



0.2 0.3 



0.4 



0.5 



1.4 
1.2 
1 



A v 0.8 



'b 0.6 
0.4 
0.2 







■5 






. .35 




Odd 




1 At- 25 





0.1 



0.2 



0.3 



0.4 



0.5 



FIG. 5. Dimensionless gap A^/Af, versus flux for different 
system sizes a) with even N , b) with odd N. 

For the case of even N, the figure (|^.a) is consistent 
with the figure 2 of ref.cl where A v /Aq was plotted in- 
stead of A v /Ab here. That figure could not show the 
interesting result that in low flux the gap increases when 
the size decreases. 

For odd N, the order parameter increases in a finite 
magnetic flux and it is maximum at tp = 1/2. This 
contradicts the arguments of ref.0 who argues that the 
suppression of the CDW order is due to pair-breaking in- 
duced by the field. There is indeed no pair-breaking effect 
since the field does not couple to the phase of the electron- 
hole pair. The field effect here is simply to change the 
position of the energy levels and thus to either reduce or 
enhance the stability of the CDW. This can be simply 
understood from the schematic Figs.((|). 

For a ring of large size, the modulation of the gap given 
by eqE^ becomes weak and harmonic. It is given by: 



A v = A 6 (l± 



1 



-2ttA 



/So 



cos 2ttp) 



which explicitely displays the exponential decrease of the 
modulation with the size of the ring. 



IV. TRANSITION TEMPERATURE 

The dependence of the critical temperature with the 
size and the flux reflects those of the gap. _ 

At the transition, A = A = 0. As in rcfE, the self- 
consistency equation can be written as, for even N: 



2 ^ 



tanh[%(n- 



¥>)] 



■ip 



Doing the same Poisson summation as above, the self 
consistency equation for the critical temperature T v is 
found to be, taking into account the parity: 

- = In 1.14^ + ]T (-l) Nm F(m/f3' v ) cos{2nmp) (11) 



m>0 



4 



where the function F is: 



F(x) = In 



cosh(27r 2 a;) + 1 
cosh(27r 2 a;) — 1 



by: 



The critical temperature of the infinite system is given 



i= In 1.14^1 
9 T b 



(12) 



so that by difference between eqs.flllD and (|12|), one has 
the flux dependence of T^/T b for different parameters 
A b = A b /2S: 

^ m>0 ^ ' 

The result is shown on Fig.fJ^). For zero flux, the vari- 
ation is very similar to that of the gap (Fig.f5j). When 
the size decreases, the critical temperature increases, in 
apparent contradiction with the figure (4) of ref.EI who 
found a decrease of the critical temperature. However 
their temperature is normalized to the gap, Ao which is 
itself size dependent, and not to What is actually 
found in refJa is an increase of the BCS ratio IS. V /T V . 
This is in agreement with our calculation of this ratio 
shctsyn in Fig.fm) and exhibits the 12% increase found in 
refl. 




FIG. 6. Variation of the critical temperature with the size, 
in zero flux. 




FIG. 7. Variation of the 

normalized BCS ratio r = (Ao/Tb)/(A{,/Tf,) with the size, 
in zero flux 

Finally, we plot in Fig.(||) the variation of the critical 
temperature with the flux, for different sizes. 
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FIG. 8. Variation of the dimensionless critical temperature 
T v /T} > with the flux, for different system sizes, when JV is even. 
The case where N is odd deduces by the same symmetry as 
in Fig. 5 

In the limit of a large system, the oscillations of the 
critical temperature become exponentially small as: 



T v = T b (l ± 2e" 1 - 14L/ « cos 2nip) 

V. PERSISTENT CURRENT 

The persistent current in the CDW phase is given by: 



y/(n + (f) 2 + A 2 



(14) 



where n G [-N/2, N/2 - 1] if N is even and n G [-(N- 
l)/2,(N- l)/2] if N is odd. 



h = — = 



25 evF 



5 



is the maximal current in one dimension. Eq.(^) for 
the persistent current is exact for a quadratic dispersion 
relatione!. The first term is the persistent current /jv(v) 
in the normal state: /jv(v) = —2Io f when N is odd and 
I N {tp) = -I (2\<p\ - 1) when N is even£3 

After summation by parts and Poisson summation, the 
flux dependence of the total current can be conveniently 
cast in the Fourier expansion which is parity dependent: 

I(cp) = 4I A V (~l) Nm Ki(2nmA ip )sin2T:rrvfi 

m>0 

K\ is a modified Bessel function of the second kindi. The 
persistent current depends on the gap which is itself flux 
dependent. For pedagogical purpose, we first show the 
flux variation of the current assuming a constant gap, 
Fig-®- 
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FIG. 9. Persistent current with a constant gap Ao, for an 
even number of particles. 

When the gap goes to zero, K\(x) p-> 1/x and one 
recovers the current of the normal statdj: 
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FIG. 10. Persistent current in the CDW phase (full line), 
for two values of the size L, in the case of even N. The 
dashed line shows the current if the gap were constant and 
the dotted-dashed line shows the current in the metallic phase. 



2 (—i) Nm 

I (if) = —Iq / sin 2Trmip 

m>0 

When the gap becomes larger than the interlevel spacing, 
the current is reduced exponentially as: 

I{cp) = (-l) Ar / e- 1 - 14i /«sin2^ 

The variation of the gap itself with the flux must be 
taken into account. Fig.(|10|a) shows the variation of the 
current with the flux for Ao = .35. When the flux in- 
creases, the gap decreases so that the current becomes 
larger (full line) than if the gap were constant (dashed 
line). When A = .25, the CDW gap vanishes at a crit- 
ical flux and the current recovers continuously its value 
in the metallic phase. This results constradicts those of 
refill who found a discontinuity in the current. When N 
is odd, the current is trivially shifted by half a period 
4>a/2, at variance with the conclusion of ref.B. 



VI. CONCLUSIONS 

We have derived the mean-field thermodynamics of a 
CDW in a small ID clean system in the presence of a 
magnetic flux. The stability of the CDW depends on the 
parity of the number N of particles. When the size de- 
creases and becomes of the order of the coherence length 
£ = fivp I V A, the CDW order parameter increases if N is 
even, it decreases if N is odd. The CDW is stabilized by 
the magnetic flux when N is odd and it is destabilized 
when N is even. These results correct those of refsfiEl 
who found that the flux always tends to suppress the 
Peierls instability. 

These are the results for a one dimensional ring. They 
can be in principle generalized to the case of a many 
channel ring. Reffl suggests that the current is simply 
multiplied by the number of chains. This is not true, as 
it is already known for the metallic phase that the cur- 
rent results from a subtle addition of the contributions 



G 



of the different channels^. Such rings with few number 
of channels can berSjpthetized using thin-film growth of 
blue bronze oxydest3'li2l. The case of a ID ring with short 
range interaction and impuritieSp-a disordered Luttinger 
liquid, has been studied recentlytil. It would be interest- 
ing to see how the discretness of the spectrum affects the 
obtained results. 

Note added in proof: After this paper was accepted, 
I have been informed by F. Von Oppen oLthe existence 
of a related work with similar conclusions!!!! Here I have 
found the analytical expressions of the harmonics expan- 
sion of the critical temperature, gap and persistent cur- 
rent.^-. They agree with previous numerical calculations. 
Rcf.ta addresses the fluctuation effect. 
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